An outburst from a massive star 40 days before a su- 
pernova explosion 
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Various lines of evidence suggest that very massive stars experience extreme 
mass-loss episodes shortly before they explode as a supernova, t 14 ' Interestingly, 
several models predict such pre-explosion outbursts. ' 5 ~ 7 1 Establishing a causal 
connection between these mass-loss episodes and the final supernova explosion 
will provide a novel way to study pre-supernova massive-star evolution. Here 
we report on observations of a remarkable mass-loss event detected 40 days 
prior to the explosion of the Type Hn supernova SN 2010mc (PTFlOtel). Our 
photometric and spectroscopic data suggest that this event is a result of an 
energetic outburst, radiating at least 6 x 10 47 erg of energy, and releasing about 



2 Ofek et al. 

10 2 M Q at typical velocities of 2000 km s 1 . We show that the temporal prox- 
imity of the mass-loss outburst and the supernova explosion implies a causal 
connection between them. Moreover, we find that the outburst luminosity 
and velocity are consistent with the predictions of the wave-driven pulsation 
model' 6 ], and disfavour alternative suggestions' 7 '. 

Type Iln supernovae (SNe Iln) are a diverse class of transient events, spectroscopically 
denned by narrow and/or intermediate- width hydrogen emission lines, up to a few thousand 
kms" 1 in width.' 8-10 ! These emission lines likely originate from optically thin circumstellar 
matter excited by the SN shock and/or radiation field. The spectra and light curves of 
these SNe are interpreted as signatures of interaction between the SN ejecta and mass that 
has been expelled prior to the explosion. ' n _13 1 The best case in which the progenitor of a 
SN Iln has been detected prior to the SN explosion' 3 ' was likely a luminous blue variable 
(LBV) - a class of objects which are known for their vigorous eruptive mass-loss events' 14 '. 
Very recently, a massive star in the nearby galaxy NGC7259, SN2009ip, had a series of 
outbursts' 15-17 ! with a typical velocity of about 500 km s _1 , possibly followed by an ongoing 
supernova event' 18-19 ! exhibiting P-Cygni lines with velocities of ~ 10 4 kms _1 . Another 
possibly related event is SN2006jc that is further discussed in SI §7. 

The Palomar Transient Factory' 20-21 ' discovered the Type Iln SN2010mc 
(PTFlOtel)' 241 in images obtained on 2010 Aug. 20.22 (UTC dates are used through- 
out this paper). It is located at a = 17 h 21 m 30. s 68, 5 = +48°07 / 47."4 (J2000.0), at redshift 
z = 0.035 which corresponds to a luminosity distance of 153 Mpc. Our ongoing search 
for precursor events in pre-explosion images of nearby PTF SNe Iln revealed a positive 
detection at the SN location in images taken prior to the SN discovery (Figure 1). Here 
we measure time relative to 2010 Aug. 20.22, which corresponds to the onset of the SN 
explosion (main event; Fig. 1). The initial bump emerged at day —37 relative to the SN 
discovery date, and peaked at an absolute magnitude of about —15 (~ 2.25 x 10 41 ergs -1 ) 
in the -Rptf band. The main SN explosion then brightened for two weeks and peaked at 
an -Rptf absolute magnitude of —18.4 (~ 5.2 x 10 42 ergs -1 ), radiating a total bolometric 
luminosity of ~ 3 x 10 49 erg, while the precursor bump radiated ~ 6 x 10 47 erg (or more 
due to the unknown bolometric correction). 
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Spectra of the supernova, showing a blue continuum with Balmer emission lines, are 

presented in Figure 2. The continuum becomes redder with time, and its slope corresponds 
to an effective temperature of over 16,000 K at day five and drops to about 8,000 K at day 
27. The Ha line has an initial width of ~ 3 x 10 3 kms" 1 at day 6, decreasing to ~ 10 3 kms" 1 
at day 14. A broad (10 4 kms _1 ) P-Cygni profile emerges by day 27. The spectra also show 
He I lines with decreasing strength, presumably due to the drop in temperature. 

The nature of the precursor bump is very intriguing and can potentially tell us a 
great deal about the SN explosion and the progenitor. The only interpretation that is 
fully consistent with the photometric and spectroscopic evidence is that the first bump 
represents an outburst from the SN progenitor about one month prior to explosion, while 
the brighter bump is initiated by a full explosion of the star a few weeks later. Below we 
analyse this model in the context of the photometric and spectroscopic data. In SI §6 we 
discuss some alternative models and conclude that they are unlikely. 

The mass ejected by the precursor burst can be estimated in various independent ways. 
By requiring that the precursor integrated bolo metric luminosity, -Eboi.prec, is lower than 
the kinetic energy of the precursor outburst (moving at velocity t> prcc ) which powers it, we 
can set a lower limit on the mass ejected in the precursor outburst M prec <; 2i? boliPrcc t>p r 2 c ps 
1.5 x 10" 2 (t> prcc /2000 kms -1 )' 2 M Q . The outburst velocity is estimated from the line widths 
of 1000-3000 kms" 1 , seen in the early-time spectra of the SN. As this mass was presumably 
ejected over a period of about one month (i.e., the outburst duration), the annual mass- loss 
rate is about 10 times higher. A similar order of magnitude argument can be used to put 
an upper limit on the mass in the precursor outburst. If some of the SN bolometric energy, 
-E'boi.sN, is due to interaction between the SN ejecta, moving at i>sn, and the precursor 
shell, and assuming high efficiency of conversion of the kinetic energy to luminosity, then 
M prcc £ 2£ boliSN ^s"N ~3x lO-VsN/lO^ms- 1 )-^. 

Another method to estimate the mass in the ejecta is based on the Ha emission-line 
luminosity and the radius of the photosphere, which is determined from a black-body fit 
to the spectra (SI Fig. 2). This method, derived in SI §2, suggests a mass-loss rate of 
<; 10" 1 M yr _1 , or a total ejected mass of ^ 10" 2 M , if we assume a month-long outburst 
with a mean wind velocity of 2000 kms" 1 . 
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Another independent method for estimating the mass in the CSM is based on the 
rise time of the supernova. Since photons diffuse through material between the SN and 
the observer, the SN rise time gives us an upper limit on the diffusion time scale and 
therefore the total intervening mass. Since the SN rise time was about one week, this 
argument suggests that if the shell is spherically symmetric its total mass cannot exceed 
about 0.4 M (see SI §5). Therefore, the kinetic-energy arguments, the Ha luminosity, the 
diffusion time scale, and the X-ray limits (see SI §3) are consistent with each other and 
indicate that the total mass lost during the outburst is of order ~ fCT 2 M . 

Our model for the sequence of events is presented in Figure 3. In a nutshell, this 
model suggests that the precursor outburst ejected ~ IO~ 2 M at a velocity of 2000 km s' 1 
about one month prior to the SN explosion. Shortly after the SN explosion, this ejected 
material was engulfed by the SN ejecta. At later times, after the optical depth decreases, 
we start seeing indications for the high velocity of the SN ejecta. We discuss some less 
likely alternative models in the SI §6. 

A surprising result is the short time between the outburst and the explosion, which 
is a tiny fraction (~ 1 0~ 8 ) of the lifetime of a massive star. Even if massive stars have 
multiple mass-loss episodes (with mass loss <; IO~ 2 M ) during their lifetime, the number 
of such episodes cannot exceed ~ 5000; otherwise, it will exceed a typical stellar mass. A 
conservative estimate shows that so far it might be possible to detect such an outburst in 
a sample of up to about 20 nearby SNe Iln, which have deep pre-explosion observations 
(see SI §7). Therefore, the probability of observing a random burst one month prior to the 
explosion is 0.1%. We conclude that such outbursts are either causally related to, or at 
least two orders of magnitude more common before, the final stellar explosion. 

When considering where such an eruption could originate, one is led to the stellar 
core, as the binding energy it can liberate by contraction and nuclear reactions is suffi- 
ciently large. Because the actual luminosity is most likely neutrino-dominated, the Kelvin- 
Helmholtz time scale for this contraction can be short enough to explain the precursor's 
rapid rise in luminosity. However, the energy liberated at the core must reach the enve- 
lope, and it therefore requires a fast and efficient mechanism' 6 ' . There are several proposed 
mechanisms that predict high mass-loss rates prior to the final stellar explosion. Recently, 
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it was suggested' 6 ] that in some massive stars the super-Eddington fusion luminosities, 
shortly prior to core collapse, can drive convective motions that in turn excite gravity 
waves that propagate toward the stellar surface. The dissipation of these waves can un- 
bind up to several solar masses of the stellar envelope. Alternatively, it was suggested' 7 ' 
that the mass loss is driven by a common-envelope phase' 25 ' due to the inspiral of a neu- 
tron star into a giant companion core (a so-called Thorne-Zytkow object), unbinding the 
companion envelope and setting up accretion onto the neutron star that collapses into a 
black hole and triggers a SN explosion shortly thereafter. This model, however, predicts 
an outflow velocity that is considerably lower than the one observed in the precursor of 
SN2010mc, and is therefore disfavoured. Another possible mechanism is the pulsational 
pair instability SN' 26 !. However, current models predict that the mass ejected in these 
events will likely be much larger than the 1O _2 M seen in SN2010mc for the observed 40 
day delay between the two last explosions. I 5 ]'' 27 ] 

The velocity and energetics of the precursor of SN2010mc are consistent with the pre- 
dictions of the wave-driven outburst mechanism' 6 ' . However, a more detailed theoretical 
investigation is required in order to test that this model is consistent with all the obser- 
vational evidence. Moreover, in SI §8, given the observed precursor luminosity, we theo- 
retically derive the mass loss to be ~ 0.05 M , in excellent agreement with our mass-loss 
estimates. Furthermore, we note that the velocity of the precursor ejecta is higher than 
predicted by other models. ^ Finally, we note that the mass-loss velocity in SN2010mc 
is considerably higher than the one observed in SN 2009ip,' 18 ' although the mass-loss rate 
was similar.' 28 ' 
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Figure 1: The light curve of SN 2010mc as obtained with the Palomar 48- 
inch telescope. The red circles are based on individual images, the squares are based 
on coadded images, while the empty triangles represent the 3a upper limits derived from 
coadded images. The error bars represent the la errors. The object magnitudes are given 
in the PTF magnitude system. I 22 !'! 23 ! Other datasets, including the Palomar 60-inch and 
Swift-XJVOT, are listed in the SI (Table 3). All the photometric and spectroscopic data 
are available via the WISeREP website' 29 '. The dashed line shows the expected luminosity 
from the radioactive decay of an ejected mass of 0.1 M of Ni 56 . Assuming that at late times 
the optical depth is sufficiently large to convert the radioactive energy to optical luminosity, 
but not too large so it will go into PdV work, this line represents an upper limit on the 
total amount of Ni 56 in the ejecta; it was set to coincide with the latest observation of the 
SN at JD 2,455,758 (see Table 3). The dotted line represents a bolometric luminosity 
equal to 50 times the Eddington luminosity for a 50 M Q star (order of magnitude estimate 
of the mass of the progenitor assuming it is a massive star). The right edges of the "S" 
symbols above the light curve indicate the epochs at which we obtained spectra (see Fig. 2). 
A full version of this light curve, including the late-time observations, is shown in the SI. 

Figure 2: Spectra of SN 2010mc, showing prominent Balmer emission lines. The 

early-time spectra are well fitted by a narrow component (~ 10 2 kms _1 ; the fit includes 
the instrumental line broadening) and a broad component extending to ~ 3 x 10 3 kms _1 at 
the first epoch, and decreasing to 10 3 kms _1 at the second epoch. At later epochs the Ha 
line develops a P-Cygni profile (similar to that of normal SNe II) with a velocity difference 
between the absorption bottom and the emission peak of ~ 10 4 kms _1 . At the first epochs 
we also detect He I lines, some of which are marked on the plot. They become weaker (see 
SI Figure 2) at later epochs, presumably due to the decrease in the effective temperature 
(see SI Figure 2). This, as well as the absence of Hell lines, indicates that our temperature 
estimate (SI Figure 2) is reasonable. 

Figure 3: Qualitative sketch of the proposed model for SN 2010mc. Panel a: 
At day ~ 0, an inner shell (purple) with a mass of ~ 1O~ 2 M , ejected about one month 
earlier, during the precursor outburst, and moving at about 2000kms _1 , is located at a 
radius of ~ 7 x 10 14 cm. An outer shell (orange), found at a large radius and moving 
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at a few hundreds kms -1 (up to 10 3 kms _1 ), was ejected at earlier times. This indicates 
that the progenitor likely had multiple mass-loss episodes in the past tens to hundreds 
of years prior to the explosion. Panel b: At day ~ 5, the SN shock front (grey line) 
moving at ~ 10 4 kms _1 is ionising the inner and outer shells which produce the broad 
and narrow Ha emission seen in the early-time spectra. Panel c: At day ~ 20, the SN 
shock engulfs the inner shell, and the intermediate width (~ 2000 km s -1 ) component of 
the Ha line disappears. Instead we detect a lOOOkms^ 1 line, presumably due to material 
ejected during previous, but likely recent, mass-loss episodes and that is found at larger 
distances from the SN. We note that inspection of the SN light curve shows that around 
day 50 there is an indication for a possible rebrightening, perhaps resulting from the SN 
ejecta colliding with such additional material ejected at earlier times. At day ~ 20, the 
photospheric temperature decreases and it becomes optically thinner, and therefore we 
begin seeing an Ha P-Cygni profile with a velocity of ~ 10 4 kms _1 . This line become even 
stronger on day 27. This reflects the unshocked ejecta below the interaction zone. 
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SUPPLEMENTARY INFORMATION 



1 Observations 

1.1 Photometry and spectroscopy 

On 2010 Aug. 25, the Palomar Transient Factory (PTF; Law et al. 2009; Rau et al. 2009) 
discovered with an autonomous machine- learning framework (Bloom et al. 2011) the Type 
Iln SN2010mc (PTFlOtel; Ofek et al. 2012c) in images obtained on 2010 Aug. 20.22 UT. 
Further inspection of the images generated by the automatic image-subtraction pipeline 
revealed that the object was detected even on Aug. 20.22 (i.e., the discovery date). 

We reanalysed the images using a careful image-subtraction pipeline. The improved 
reduction clearly shows a variable source at the location of the SN prior to the discovery 
date. Table 1 presents the photometric measurements of SN 2010mc, including observations 
made with the Oschin Schmidt 48-inch telescope, the Palomar 60-inch telescope (Cenko 
et al. 2006), and the Swift-XJVOT (Gehrels et al. 2004). We note that the ground-based 
observations were processed using image subtraction. 

A version of the light curve showing the late-time observations of the SN is presented 
in SI Figure 4. 
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Figure 4. The light curve of SN2010mc (PTF lOtel) as obtained with the Palomar 48-inch 
telescope, including late-time observations. See main text Figure 1 for details. 
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Table 1: SN2010mc (PTF lOtel) photometric mea- 



surements. Time is measured relative to JD = 
2,455,428.722 day. A_ and A + are the lower and up- 
per time range in which the photometry was obtained. 
These columns are indicated in cases where we coadded 
data from several images. N is the number of images. 
PTF magnitudes are given in the PTF magnitude system 
(Ofek et al. 2012a; Ofek et al. 2012b). The Swift/XJVOT 
measurements represents the aperture magnitude at the 
position of the SN and include the host-galaxy light. The 
two last UVOT measurements were taken at late times; 
thus, they can be regarded as the magnitudes of the host 
galaxy. Based on them, we conclude that in order to 
correct the early-time measurements for the host-galaxy 
contamination, about 0.05 mag should be added to the 
SN magnitudes. 
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UT Date 


Telescope 


Instrument 


Ha flux 


broad Ha vel. 


narrow Ha vel. 








(erg cm~ 2 s _1 ) 


(kms- 1 ) 


(kms -1 ) 


2010 Aug. 26 


Gemini 


GMOS 


9.2 x 10~ 15 


3400 


800 


2010 Sep. 03 


Lick 


Kast 


1.2 x 10~ 14 


1100 


160 


2010 Sep. 10 


KPNO 


RC spectrograph 


9.2 x 10~ 15 


670 


110 


2010 Sep. 16 


Lick 


Kast 


1.4 x 10~ 14 


4000 


230 



Table 2. Log of spectroscopic observations of SN2010mc (PTF lOtel) and properties of the Ho 
emission line. The fluxes and velocities were measured by fitting two-component Gaussians which 
were convolved with the instrumental broadening as estimated from atmospheric emission lines. 



In addition, we obtained spectra of SN2010mc on four epochs, as listed in Table 2 
and presented in Figure 2 of the main paper and SI Figure 5. In order to measure the 
line fluxes in a consistent way, we first calibrated the flux of each spectrum such that the 
synthetic -Rptf magnitude derived from the spectrum was identical to the -Rptf magnitude, 
interpolated to the epoch of each spectrum, from the PTF magnitudes listed in Table 1. 
All of the spectra are available online from the WISeREP website* (Yaron & Gal- Yam 
2012). Each spectrum was smoothed using a median filter and a black-body curve was 
fitted. The derived effective temperatures and radii are presented in SI Figure 6. 

The X-ray observations of PTF lOtel are presented and discussed in Section 3. 

2 Derivation of mass loss as a function of the Ha line luminosity 

Here we derive a relation between the mass in the circumstellar matter (CSM), the Ha 
emission-line luminosity Ln a , and the radius r. We assume that the CSM has a wind- 
density profile of the form p = Kr~ 2 , where K = M / (4nv w ) , M is the mass-loss rate, and 
v w is the wind/outburst velocity. In this case the particle density profile is given by (e.g., 
Ofek et al. 2012c) 

1 M 
n ~ 

(/i p ) 4irm p v w r 2 

* 3.02X10-JL ^— -( v - r rY^— r 2 cm-3, (1) 

(^O.lMoyr" 1 VlOkms -1 / Vl0 15 cm/ ' 1 ; 

*http://www. weizmann.ac.il/astrophysics/wiserep/; Weizmann Interactive Supernova (data) 
REPository. 
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where (/x p ) is the mean number of nucleons per particle (mean molecular weight). For our 
order-of-magnitude calculation, we adopt (/i p ) = 0.6. Assuming the SN radiation field can 
ionise all of the hydrogen in the CSM, the mass of the hydrogen generating the Ha line is 



M H ~ h off • ( 2 ) 



Here vu a is the Ha frequency (4.57 x 10 14 Hz), a^ a is the Ha effective recombination 
coefficient, ~ 8.7 x 10" 14 T~° 89 cm 3 s" 1 (Osterbrock & Ferland 2006), where T 10k is the 
temperature in units of 10,000 K. We assume T 10k = 1. In a wind profile, the integrated 
mass from radius r to r is 

M = j Anr 2 Kr~ 2 dr = ArcK(r - r ) « AnKr, (3) 

where the last step assumes that r <r. By substituting Equation 1 into Equation 2 and 

setting it equal to Equation 3, we can get a relation between the mass-loading parameter 

K, the Ha luminosity, and the radius: 

/ (// p )mgL Ha r \i/2 
^ V 47rhua cG J 

~ 7.1xl0 15 f ^ Ha - Y /2 (—^- V^gcm- 1 . (4) 
VlO^ergs- 1 / Vl0 15 cm/ s v ; 

The reason for the inequality is that in practice, it is possible that not all of the hydrogen is 
ionised or that the temperature of the gas is too high (i.e., a^ depends on temperature). 

Assuming the outburst ejected material with a velocity of 2000 km s" 1 at day —37, 
then at the time the first spectrum was taken (day 6), the ejecta radius is r w 7 x 10 14 cm. 
This radius is consistent with the one derived from the black-body fit to the spectrum at 
the first epoch (see SI Figure 6). Using the measured flux of the Ha broad component 
at day 6 (9 x 10~ 15 ergcm~ 2 s^ 1 ; Table 2), the luminosity of the Ha broad component is 
~ 3 x 10 40 ergs -1 . Substituting these values into Equation 4, and assuming a wind velocity 
of 2000 km s -1 , we find a mass-loss rate of > 10~ 1 M Q yr -1 . Assuming a one month long 
outburst this gives a total mass of about 10~ 2 M Q . This value is in excellent agreement 
with the mass-loss rate we independently find based on the kinetic-energy arguments and 
the diffusion time scale. 

We note that this derivation assumes that the narrow lines are due to the interaction 
of the CSM with the SN radiation field. One caveat, is that it is possible that some, or all, 
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UT Date 


Time since disc. 


Exp. Time 


2<T upper limit 




(day) 


(s) 


(ctks- 1 ) 


2010 Aug. 27.0 


6.8 


1342 


2.2 


2010 Aug. 31.4 


11.2 


2256 


1.3 


2012 July 11.7 


691.5 


4818 


0.62 



Table 3. Swift-XRT observations of SN2010mc (PTFlOtel). 

of the line flux originates from a shock interaction or from a slow moving ejecta. However, 
the order of magnitude consistency between this mass-loss estimate and other mass-loss 
estimators suggest that such an effect is not significant. 

3 X-ray observations 

The Swift-XRT observations of SN2010mc (PTFlOtel) in the 0.2-10 keV band are listed 
in Table 3. For each Swift-XRT image of the SN, we extracted the number of X-ray counts 
in the 0.2-10 keV band within an aperture of 7.2" (3 pixels) radius centred on the SN 
position. We chose a small aperture in order to minimise any host-galaxy contamination. 
We note that this aperture contains ~ 37% of the source flux (Moretti et al. 2004). The 
background count rates were estimated in annuli around each SN, with an inner (outer) 
radius of 50" (100"). Assuming a photon spectrum n p h oc E~ r with T = 2, a Galactic 
neutral hydrogen column density of 2.34 x 10 20 cm -2 in the direction of the SN (Dickey 
& Lockman 1990), and correcting for the photometric aperture losses, we find a 2cr upper 
limit on the luminosity in the second epoch of 3.5 x 10 41 ergs -1 . 

Following Immler et al. (2008), the X-ray emission from the optically thin region is 
given by 

/oo 
Anr 2 A(T)n 2 dr, (5) 

where A(T) is the effective cooling function in the 0.2-10 keV range. Assuming an optically 
thin thermal plasma with a temperature in the range 10 6 -10 8 K (Raymond et al. 1976), 
we adopt a value of A(T) « 3 x 10~ 23 ergcm 3 s -1 . Substituting Equation 1 into Equation 5 

we get 

L x ~ 47rA(T)^ 2 ^ r 
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~ 2.4xl0 40 f ,W ^ rV-^—rV^Jergs- 1 . (6) 

VO.OlM yr- 1 / V 2000 km s -1 / Vl0 15 cm/ & V ; 

Here, we added a correction factor due to absorption in the wind, where r and r bf are the 
Thomson and bound-free optical depths, respectively. In our case the Thomson optical 
depth will be 1 for M w O.lMoyr" 1 (e.g., Ofek et al. 2012c). Although the Thomson 
optical depth is well known, when the optical depth is of the order of a few, Compton scat- 
tering is expected to reprocess more energetic photons into the 0.2-10 keV band (Chevalier 
& Irwin 2012; Svirski et al. 2012). Since the exact X-ray spectrum is not known (Katz et 
al. 2011; Svirski et al. 2012), a proper calculation of L x when r 1 is not straightfor- 
ward. Regarding the bound-free optical depth, although it is likely that all the hydrogen 
is ionised, Chevalier & Irwin (2012) argued that if the shock velocity is below ~ 10 4 kms _1 
some of the metals will be neutral. In this case bound-free absorption will be important 
and it is given by (Ofek et al. 2012c) 

Tbf ~ VO.OlMoyr" 1 / V2000kms _ V U0 15 cm/ VlkeW ' U 

where Ex is the X-ray energy. 

If we ignore the Thomson optical depth (i.e., r = 0), and if we assume a flat spec- 
trum (i.e., constant number of photons per unit energy) and the Swift-XBT spectral re- 
sponse, from Equations 6 and 7 we find that the mass-loss rate is either smaller than 
about 0.1M Q yr~ 1 or larger than about lM yr _1 . The low bound of M < O.lMQyr" 1 is 
consistent with the other mass-loss estimates. 



4 Radio observations 

We observed SN2010mc (PTF lOtel) with the Jansky Very Large Array (JVLA f ) on 2012 

Sep. 10, about two years after the SN discovery. The observation was undertaken in 

the K band (21.8 GHz). We used J1727+4530 for phase calibration and 3C286 for flux 

calibration. The data were reduced using the AIPS software. We did not detect SN 2010mc 

with an upper limit of 0.2 mJy (3a). 

^The Jansky Very Large Array is operated by the National Radio Astronomy Observatory (NRAO), a 
facility of the National Science Foundation operated under cooperative agreement by Associated Univer- 
sities, Inc. 
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5 Diffusion time scale limit on mass loss 



Another observable that can be used to constrain the mass in the CSM is the rise time of 
the SN light curve. If a considerable amount of material is present between the SN and 
the observer, then photon diffusion will slow down the rise time of the SN light curve. 
Therefore, the maximum observed SN rise time can be used to put an upper limit on the 
amount of mass between the SN and the observer. The diffusion time scale in a wind 
profile is given by (e.g., Ginzburg et al. 2012) 
nK ^ ( c \ 

Here i> e j is the SN shock velocity. In our case, assuming the rise time of the SN is about 
seven days, we can use Equation 8 to set an upper limit, M < 0.4M Q yr -1 . 

6 Alternative interpretations 

Here we consider alternative interpretations regarding the nature of the precursor bump 
and show that they are not consistent with the observations. One possibility is that this 
feature is due to a shock breakout (e.g., Colgate 1974) prior to the SN rise. However, 
given the feature's luminosity, its duration is at least an order of magnitude too long to be 
consistent with a shock breakout (e.g., Nakar & Sari 2010; Rabinak & Waxman 2011). 

A second possibility is that the first bump represents a SN explosion while the second 
peak is due to the interaction of the SN ejecta with dense CSM and the conversion of the 
SN kinetic energy into radiated luminosity. However, this scenario can be rejected based 
on the fact that by fitting a black-body spectrum to the spectra (see SI Figure 6), we can 
find the temperature and photospheric radius evolution between days 5 and 27. These fits 
suggest that the photosphere is expanding with a mean velocity of >4300kms _1 , and that 
the black-body radius at day 6 is ~ 6 x 10 14 cm. Extrapolating backward, this suggests 
that the SN outburst happened after day —10. 

Another argument against this scenario is that the SN spectrum at day 27 after discov- 
ery (main text Figure 2), shows a P-Cygni profile with a velocity of ~ 10 4 kms -1 . At the 
same time, the black-body radius of the SN photosphere is ~ 1.5 x 10 15 cm. These values 
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are roughly consistent (to within 30%) with an explosion time at day 0, and are less con- 
sistent with explosion at day —37 day relative to discovery. A caveat in this statement is 
that the photospheric radius which controls the black-body temperature of the continuum 
is not necessarily at the same position where the P-Cygni absorption is generated. 

7 The Type Ibn SN 2006jc 

Another possibly related event is SN2006jc (Nakano et al. 2006). Unlike SN2010mc 
(PTF lOtel) and SN 2009ip, this was a hydrogen-stripped Type Ibn SN. In this case an 
optical outburst was observed in 2004, about two years prior to the SN explosion (Nakano 
et al. 2006), creating a dense circumstellar shell (Foley et al. 2007; Immler et al. 2008; 
Smith et al. 2008; Ofek et al. 2012c). Foley et al. (2007) suggested that the progenitor 
was a Wolf-Rayet star that had an outburst two years prior to its explosion, creating the 
CSM; it may have recently transitioned from the LBV phase. However, Pastorello et al. 
(2007) argued that SN2006jc took place in a binary system in which one star (a Wolf- 
Rayet star) was responsible for the SN explosion while the companion was a luminous blue 
variable (LBV) that generated the eruption in 2004. Our argument regarding the causality 
between pre-explosion outbursts and the SN explosions disfavours the explanation given by 
Pastorello et al. (2007) for SN2006jc. However, for SN2006jc, the time difference between 
the precursor event and the explosion is about two years (rather than one month), rendering 
the argument a factor ~ 20 weaker than in the case of PTF lOtel. 

8 The probability of detecting a precursor outburst prior to the explosion 

The specific precursor event reported in this paper was found as a part of a search for 
outbursts among type-IIn SN in PTF data. The sample we used included only four nearby 
type-IIn SNe with good pre-explosion coverage. Nevertheless, in estimating the probability 
of detecting a precursor prior to the explosion we assumed a sample size of 20 (five times 
higher than our actual sample size), accounting for other samples in which outburst events 
could have been detected (see details below). This sample size thus facilitates a conservative 
rate estimate. 

Thus far the brightest precursor outbursts (e.g., PTF lOtel) have an absolute magnitude 
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of about —16. Using previous-generation transient surveys, such bursts can be detected 
up to a distance of ~ 100 Mpc, while the new-generation surveys (e.g., PTF) may be able 
to detect them to somewhat larger distances. Among the 165 SNe Iln listed in the IAUC 
website, about 70 have distances below 100 Mpc. Many of these SNe were found using small 
telescopes which are limited by their depth. Based on the PTF experience, we estimate 
that not more than 20% of these SNe have sufficiently good pre-explosion observations that 
will allow the detection of a precursor prior to the explosion. In addition, PTF has already 
found 91 SNe Iln (some listed in the IAUC website), of which only a few are sufficiently 
nearby, and have high cadence and almost continuous observations prior to the explosion. 
Thus, we conclude that up to 20 SNe Iln discovered so far have good enough pre-explosion 
observations to find precursor outbursts. 

We note that even if we assume that all of the ~ 250 SNe Iln discovered so far have 
sufficiently good observations prior to the SN explosion, this will not change our conclusion 
that precursor outbursts are more common, or are causally connected with the final stages 
of massive-star evolution. 

Finally, we note that we cannot completely rule out the possibility that the SN and 
precursor are two unrelated events (e.g., explosions of two different stars in a star cluster). 
However, we think this is an unlikely possibility given the evidences presented in this 
paper, including consistency of the mass-loss estimators with the properties of the outburst 
(luminosity and time before SN), and the small number of type-IIn SNe in which we 
searched for precursor events. 

9 Theoretical determination of the precursor mass loss 

The precursor is characterised by a rapid brightening of the radiative flux which is sustained 
at a highly super-Eddington level for a month. The relatively short time scales imply that 
some efficient mechanism exists to very quickly transport the flux out, even faster than 
convection can. The mechanism suggested by Quataert & Shiode (2012), where the energy 
is transported outward by acoustic waves, offers a natural explanation. However, this 
mechanism can operate only in regions with a sufficiently high density. At a region where 
the density, p, satisfies L pa L convJnax pa 4nr 2 pv^ (where v s is the speed of sound), the 
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acoustic energy flux L is converted into a radiative flux L. The mass-loss rate expected 
in such a case is M 47rr 2 pt> s ^/fg- Assuming t> s ps 60kms _1 at the base of the wind 
gives roughly 8 M for the mass ejected during the precursor. A similar problem exists for 
the great (super-Eddington) eruption of rj Carina. The naive expectation for the mass loss 
is orders of magnitude larger than what was actually observed (Shaviv 2000). 

To overcome this discrepancy, it was suggested that instabilities reduce the effective 
opacity of super-Eddington atmospheres. This implies that a wind is accelerated only from 
lower-density regions where a typical scale height in the atmosphere becomes of order unity 
and the effective opacity returns to its microscopic value. Under these assumptions, the 
mass-loss was determined by Shaviv (2001) to be AM = WAT(L — L Edd )/(cv s ), where AT 
is the duration of the outburst and L is the luminosity at the base of the wind (where part 
of the latter is converted into mechanical energy of the mass loss). W is a dimensionless 
constant which depends on the geometry of the instabilities, and empirically was found to 
be around 3 (±0.5 dex). 

We can use this expression to predict the expected mass loss from the super-Eddington 
precursor. To do so we have to add the mechanical energy of the wind to the observed 
luminosity and obtain the total luminosity at the base of the wind. The mechanical energy 
itself is given by L mec h = AM/2(v% + v^ sc ). Here is the terminal velocity at infinity 
and f csc is the escape velocity. We crudely take t> esc ~ v^. Since the luminosity at the 
base is significantly super-Eddington, we can ignore the L Edd term. If we further consider 
that v s ~ 30kms _1 , we obtain AM ps 0.05 M q (with an uncertainty of about 0.5 dex). 
This value sits comfortably in the observationally determined range of possible mass loss 
in SN2010mc (PTFlOtel). 
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Figure 5. Close-up view of the Ha line in the spectra of SN2010mc (PTFlOtel). The best- 
fit width and total flux of a two-component model at each epoch are listed in Table 2. The 
width of the Ha emission is decreasing significantly after the first epoch, from 3400 km s -1 to 
^ 1000 km s -1 . In the last epoch a P-Cygni profile develops, with a velocity difference between 
the peak of the emission and the bottom of the absorption of over 10 kms . The velocity of 
the Ha absorption component in the P-Cygni feature is marked with a 12,000 kms -1 label. The 
P-Cygni profile is better seen in the stretch of Figure 2 (main text). Also visible is the decrease 
in strength of the Hel lines, presumably due to the decrease in temperature (see SI Figure 6). 
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Figure 6. Black-body temperature and radius derived from the spectra of SN2010mc 
(PTFlOtel). The temperature estimate should be regarded as a lower limit, because the peak 
of the emission is outside the optical region and metal line blanketing can reduce the flux in the 
blue parts of the spectrum. The reduction in temperature is consistent with the fact that the He I 
lines are becoming weaker at later epochs (SI Figure 5). We note that the best-fit temperature at 
day 6 (near the first spectroscopic observation) based on the broad-band Swift/WOT and P60 
photometry is about 19,000 K. The nondetection of He II lines also suggests that the temperature 
is not much higher than the value derived here. 
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